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Abstract 



The fragmentation resulting from peripheral Au + Au collisions at an 
incident energy of E = 35 MeV/nucleon is investigated. A power-law charge 
distribution, A~'^ with r 2.2, and an intermittency signal are observed 
for events selected in the region of the Campi scatter plot where "critical" 
behavior is expected. 



PACS : 05.70.Jk, 25.70.Pq, 64.70.Fx 
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Initiated by the observation of fragments in the final stages of the reaction exhibiting a 
power law in fragment charge distributions [|I| , and stimulated by the similarity of the nuclear 
matter equation of state with that of a van der Waals gas [0 , the possibility of observing a 
liquid-gas phase transition in nuclear systems has been the subject of intensive investigations 
This interest increased recently with attempts of extracting critical exponents of 
fragmenting nuclear systems produced in Au + C collisions at E = 1 GeV/nucleon ||5| 
and a "caloric curve" for projectile fragmentation reactions in Au + Au collisions at 600 
MeV/nucleon [0. In this letter, we report results obtained for peripheral Au + Au collisions 
at E = 35 MeV/nucleon which display some characteristics similar to those predicted for 
near-critical systems. 

The experiment was performed at the National Superconducting Cyclotron Laboratory 
of the Michigan State University. Fragments with charge up to Z = 83 were detected at 
3" < Oiab < 23° by the Multics array [H, and charged particles with charge up to Z = 20 
were detected at 23° < 9iab < 160° by 159 phoswich detector elements of the MSU Miniball 
1^. The charge identification thresholds were about 1.5 MeV/nucleon in the Multics array, 
independently of the fragment charge, and about 2, 3, 4 MeV/nucleon in the Miniball for Z 
= 3, 10, 18, respectively. The geometric acceptance of the combined apparata was greater 
than 87% of An. 

Our analysis is guided by calculations P,|lO|,|ll| which predict that the "critical" excitation 



energy may decrease when the system is either charged and/or rotating. Short- lived systems 
formed in central Au + Au collisions are predicted to expand and undergo a multi-fragment 
breakup due to the high combined charge of projectile and target nuclei. Evidence for 



such Coulomb driven breakup of a single source has, indeed, been observed |]T2[ in central 
collisions for this reaction. For larger impact parameters, however, several smaller sources 
can emerge corresponding to the decay of projectile and target-like residues and a "neck" |T^ 



that momentarily joins them. With an appropriate reaction filter, one might then hope to 
select primary fragments with excitation energies. Coulomb charges, and angular momenta 
appropriate to bring the system into different portions of the instability region |1T0| , |l4 |. 



Studies with the Classical Molecular Dynamics (CMD) model indicated that critical 
behavior may be achieved for peripheral Au + Au collisions at E = 35 MeV/nucleon ||15|| . 
When these CMD results are filtered by the acceptance of the Multics-Miniball arrays, the 
signals which may be indicative of criticality become washed out due to poorly detected 
events, but they can be recovered by restricting the analysis to more completely detected 
peripheral events for which the largest projectile-like fragment (PLF) is detected. Guided 
by these calculations, we select events for which the largest fragment has a velocity along 
the beam axis larger than 75% of the beam velocity and the total detected charge is between 
70 and 90. For these events, the total detected linear momentum is larger than 50% of the 
beam momentum. Their distribution of charged particle multiplicities, N^, is shown by the 
solid histogram in Fig. 1. 

Figure 2 shows a scatter plot of ln{Z^^^^) versus ln{M2) ("Campi scatter plot" [p!6| ) 
where .^4ax is the charge of the heaviest fragment and Ml is the second conditional moment 
of the charge distribution detected in the j-th event. 



M!f^ = Y.Z%{Z) (1) 
z 



Here, nj{Z) denotes the number of fragments of charge Z detected in the j-th event, and the 
summation is over all fragments but the heaviest detected one. Theoretical investigations 
suggest that such plots may be useful in characterizing near-critical behavior of finite systems 
T6| . The calculated Campi scatter plots typically exhibit two branches: an upper branch 



with a negative slope containing largely undercritical events (e.g. T < Tcru in a liquid-gas 
phase transition or p > pcnt in a percolation phase transition) and a lower branch with 
a positive slope containing largely overcritical events (T > Tent or p < Pcrit)- The two 



branches were shown to meet close to the critical point of the phase transition [|6|, p!6| , p!7 . 

The data shown in Fig. 2, display two branches similar to the ones predicted for under- 
critical and overcritical events. In the top-right part, close to the intersection of these two 
branches (potentially containing near-critical events), a separate island is observed which is 
due to fission events, as first noted by Ref. Jl^. By appropriate gates in the Campi plot. 



these fission events are removed from tlie following analysis. 

To further investigate the two branches observed in Fig. 2 and the region where they 
intersect, we employ three cuts selecting the upper branch (cut 1), the lower branch (cut 3) 
and the intersection region (cut 2); these cuts are indicated in Fig. 2. The charged particle 
multiplicity distributions observed for these three cuts are shown as dashed histograms in 
Fig. 1. Cuts 1 and 3 largely select low and high multiplicity events corresponding to very 
peripheral and central collisions (assuming on the average a monotonic relation between 
Nc and impact parameter); cut 2 represents a wide range of charged particle multiplicities 
and thus may involve a wide range of intermediate impact parameters. Thus emission from 
a unique source cannot be ascertained for cut 2, and it is likely that this cut contains 
contributions from projectile and target-like sources and from the neck [|13| which emits 
lighter fragments {Z ?5i 6 — 9) with enhanced probability as compared to the projectile residue 
1131 . However, one does not exclude that this large multiphcity distribution is related to the 
occurrence of large fluctuations as expected at the critical point. 

Fragment charge distributions, not corrected for detection efficiency, are presented [^] for 
the three cuts in Fig. 3. Cut 1 (crosses) contains both light fragments and heavy residues and 
thus resembles the "U" -shaped distributions predicted by percolation calculations in the sub- 
critical region ||18|. For cut 3 (open circles), one observes an unusually fiat charge distribution 
similar to the one previously reported |jl2| for central collisions which were selected without 
the specific constraints employed in this paper and attributed to a Coulomb-driven breakup 
of a very heavy composite system ||12| (The steep fall-off at large Z is an artifact of the 
selection of events with Z = 70 — 90 used in this paper). For cut 2 (solid points), a fragment 
charge distribution is observed which resembles a power-law distribution, P{Z) oc Z^'^, with 
r ^ 2.2. For macroscopic systems exhibiting a liquid-gas phase transition, such a power-law 
distribution is predicted to occur near the critical point ||2^ . However, it is not yet known by 



how much the final fragment distributions differ from the primary ones after the sequential 
decays of particle unstable primary fragments. 

Figure 4 shows [jl9| the logarithm of the scaled factorial moments (SFM), defined as 



^*(^^) = ^z,../ss , (2) 



{i = 2, 5), as a function of the logarithm of the bin size Ss. In the above definition 
of the SFM, Ztot = 158 and i is the order of the moment. The total interval [1, Zfot] is 
divided into M = Ztot/Ss bins of size 6s, rtk is the number of particles in the fc-th bin 
for an event, and the brackets <> denote the average over many events. A linear rise of 
the logarithm of the SFM versus 5s (i.e. Fi oc 5s~^^) indicates an intermittent pattern of 
fluctuations [|T^pT],p2|. Even though this quantity is ill defined for fragment distributions 



several theoretical studies have indicated that critical events give a power law for 
the SFM versus the bin size []6|,|2^. For cut 3 (right part of the figure), the logarithm 
of the scaled factorial moments ln{Fi) is always negative (i.e. the variances are smaller 



than Poissonian [^) and almost independent on 5s] there is no intermittency signal. The 
situation is different for cut 2 (central part). The logarithm of the scaled factorial moments 
are positive and almost linearly increasing as a function of —ln{5s), and an intermittency 
signal is observed. Region 1, corresponding to evaporation, gives zero slope. Increasing or 
reducing the size of the three cuts in the respective regions does not change significantly these 
results. The interpretation of experimentally observed intermittency signals may, however, 
be problematic due to ensemble averaging effects |Q, even though calculations show that 



impact parameter averaging only increases the absolute value of the SFM [^. Since cut 2 
may involve a large range of impact parameters, the observed intermittency signal could be 
an artifact of ensemble averaging and can, therefore, not be taken as a definitive proof of 
unusually large fluctuations in a sharply defined class of events. 

In conclusion, we have analyzed fragment production in Au + Au collisions at E = 35 
MeV/nucleon. Events were selected by requiring a total detected charge between 70 and 90 
and the velocity of the largest detected fragment larger than 75% of the projectile velocity. 
A Campi scatter plot of these events displays two branches similar to the sub- and overcrit- 
ical branches observed in theoretical studies. The selection of events from the intersection 
of these two branches (which has been associated with critical events in theoretical studies) 
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shows a power law charge distributions with an exponent of r ~ 2.2 similar to that charac- 
terizing the mass distribution near the critical point of a liquid-gas transition. These events, 
further, display an intermittent behavior similar to that expected for near-critical events. 
While these signatures have been associated with near-critical events, we must caution that 
the effects of finite experimental acceptance and event mixing with possible contributions 
from the decay of projectile-like fragments and the neck- region are not yet sufficiently well 
understood to allow an unambiguous conclusion of critical behavior in the present reac- 
tion. Our work does, however, show that different regions of the nuclear phase diagram can 
be probed at one incident beam energy by selecting events according to different impact 
parameters and/or energy depositions. 
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FIGURES 

FIG. 1. Multiplicity distribution of the events selected for the analysis (solid histogram). The 

three dashed histograms (1, 2 and 3) represent the multiplicity distributions of the events falling 
in the three cuts drawn in Fig. 2. 

FIG. 2. Campi scatter plot, ln{Zjnax) versus ln{M2). The three different regions are discussed 
in the text. Fission events are to the right of region 2. 

FIG. 3. Charge distributions for the three cuts indicated in Fig. 2. The curve represents a 
power law distribution Z~'^ with r = 2.2. 

FIG. 4. Scaled factorial moments ln{Fi) versus —ln{5s) for the three cuts indicated in Fig. 2; 
the left part of the figure corresponds to cut 1, the central part to cut 2 and the right part to cut 
3. Solid points represent the SFM of order i = 2, open circles i = 3, open squares i = 4, and open 
triangles z = 5. 
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